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ABSTRACT 

Based on an analysis of data obtained with the Wide Field Camera 3 (WFC3) 
on the Hubble Space Telescope (HST) we report the identification of two distinct 
stellar populations in the core of the giant Hn region 30 Doradus in the Large 
Magellanic Cloud. The most compact and richest component coincides with the 
center of R136 and is ~ 1 Myr younger than a second more diffuse clump, located 
~ 5.4 pc toward the northeast. We note that published spectral types of massive 
stars in these two clumps lend support to the proposed age difference. The mor- 
phology and age difference between the two sub-clusters suggests that an ongoing 
merger may occurring within the core of 30 Doradus. This finding is consistent 
with the predictions of models of hierarchical fragmentation of turbulent giant 
molecular clouds, according to which star clusters would be the final products of 
merging smaller sub-structures. 

Subject headings: Magellanic Clouds — stars: imaging — galaxies: star clusters: 
individual — galaxies: star formation 
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Introduction 



As the most lum inous and best known massive star-forming region in the Local Group 
(LG. lKennicuttlll984[ ). the Tarantula Nebula in the Large Magellanic Cloud (LMC) offers us 
the unique opportunity to investigate the process of star formation (SF) in an environment 
that, in many ways (such as metallicity, dust content, and star formation rate), resembles 
the extreme conditions of the early universe and distant star-forming regions. The star 
formation history ( S FH) o f Tarantula's ionizing cluster NGC 2070 (a.k.a. 30 Dor) is complex. 
Walborn fc Blades! ( 119971 ) defined 30 Dor as a two-stage starburst, in which the radiation 
from the compact core R136 has triggered a new generation of stars at a projected distance 
of ~ 1' (corresponding to ~ 15 pc at the the distance of the LMC). 

The complexity of NGC 2070 is reflected also in its atypical lumin osity profile, which 
shows a break at ~ 10" (~ 2.5 pc) and a bump at ~ 30" (~ 7pc; IMackey &: Gilmore 



20031 ). More recent observations in the near IR with Multiconjugate Adaptive Optics (MAD) 
indicate that the innermost ~ 30" of the luminosity profile can be fitted with a single 



Elson. Fall, fc Freeman! (119871. EFF) pow erlaw, once the asymmetric shape of 30 Dor is 
taken into account ( jCampbell et al.ll2010l ). 



Here we take advantage of the high spatial resolution and sensitivity of the Wide Field 
Camera 3 (WFC3) on board of the Hubble Space Telescope (HST) to study the innermost 
~ 2'44" x 2'44" of NGC 2070. In particular, we use public UVIS and IR channel observa- 
tions (P.I. R.W. O'Connell, GO-11360) to characterize the spatial distribution of its stellar 
populations. 



Data 



Deep WFC3 observations of R136, the core of NGC 2070, were retrieved from the 
Mikulski Archive for Space Telescopes (MAST) . The data-set consists of 6 broadband filters 
(F336W, F438W, F555W, F814W, F1 10W, and F160W), ran ging from the near-UV to the 
near-IR, and is described in detail in |Pe Marchi et al.l (1201 lh . The images were processed 
through the standard calibration pipeline CALWF3, using the most up-to-date calibration 
files. 

Positions and fluxes of the stars in the WFC3/UVIS images were measured using 
img2xym_wfc3uv, while we used img2xym_wfc3ir to measure posit ions and fluxes in th e 
WFC3/IR images. Both codes were adapted from img2xym_WFI (jAnderson et al.ll2006l ). 
The astrometry and photometry of WFC3/UVIS da t a wer e cor rected for pixel area varia- 
tions and geometric distortion as in iBellini fc Bedinl (120091 ) and iBellini et al.l ( 1201 ll ). while 
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for the WFC3/IR data we used the solutions from Anderson et al. (in preparation). Stars 
were detected independently in each filter. The WFC3 photometric catalogs were calibrated 
into the Vegamag photometric system using the most up-to-date values of the zeropoints, 
available at the WFC3 webpagejj, and then combined into a single catalog. Approximately 
3000 stars with photometric errors smaller than < 0.15 mag were found in all six bands. 



3. Analysis 



Fig. [TJ shows three color-magnitude diagrams (CMDs) of the stars found in the region 
covered by the UVIS observations for different combinations of the filters F336W, F438W, 
F555W, and F814W. In each plot, stars were selected to h aye photometric erro rs < 0.15 
in both the plotted filters. Padova isochrones (dashed lines; iMarigo et al. Il2008|) for main 



seque nce (MS) and evolved stars, and Pisa isochrones (continuous lines; iTognelli et al. 
201 lh for pre-main sequence (PMS) sta rs, are superimposed to guide the eye. In doing 
this we assumed a distance of 50 kpc ( jSchaefer et al.l 120081 ). and a metallicity Z=0.008 



( IHelling. Winters. &: Sedlmayrl 120001 ). Isochrones are also corrected for reddening. 



Inspection of these CMDs shows that the stellar content is complex: the bright and blue 
upper-MS (i.e. mFsu < 21; mF55sw — "T,F8i4 J$ 0.6 in panel C) is composed of intermediate 
and high mass stars, but at redder colors (i.e. mpsssw — ^fsuw ^ 1-5 in panel C) the 
majority of the sources are PMS stars in the mass range 0.5 < M < 5.0 M and ages 
> 1 Myr. Pisa isochrones show that PMS stars become fainter and bluer with age, to the 
point that after ~ 10 Myr they are indistinguishable from the low-mass MS (mpsMW — 21; 
^F555W - m F8U < 1.0 in panel C). 

Fig. [2] shows the stellar density distributions in the ~ 2'44" x 2'44" area around R136 in 
the filters F336W, F438W, F555W, F184W, and F160W. In each map the pixel scale is 2", 
and the density levels are the square root of the number of stars. The spatial distribution of 
the stars in the F336W, F438W, and F555W maps is remarkably similar, with the majority 
of the stars "clustered" at R.A. = 05 h 38 m 42.6 s , Dec = -69°06'03", J2000, which coincides 
with the center of 30 Dor (Figs. [3] and HJ). The apparent decrease of star counts in the 
F814W and F160W maps where R136 is located is caused by the high crowding, since no 
completeness corrections were applied to the maps. 

In each map the density contours are not spherical. On the contrary they show an arc 
that departs from R136 to the Northeast. This feature is clearly visible in all the filters, from 



http: / /www. stsci.edu /hst /wfc3 /phot_zp_lbn 
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the UV to the near IR, making it unlikely that local high reddening can be at the origin of 
this structure. The low impact of variable reddening in the stellar counts is also shown by 
the ratio between the F160W and F555W isodensity maps (Fig. [2]- lower-right panel). 

In Fig. [5] the (mpsssw, ^F555W — ^F8i4w) CMD of the stars found within 10" from the 
center of R136 is compared to that of the stars found in a 14.4" x 27.4" box centered on the 
Northeast clump (see also Fig.H])s. In both plots, most of the brighter stars (mp555w ^ 19.5) 
trace a well defined blue (mp555w — m F8i4W IS 0-5) upper-MS, while below this value the 
majority of the stars show the red colors typical of intermediate and low-mass (0.5 < M < 
5.3 Mq) PMS stars. Although the Northeast Clump is less then ~ 22" (corresponding to 
a projected distance of ~ 5.4 pc) away from R136, the redder and broader upper-MS of 
R136 indicates that the two populations are affected by different dust extinction. Isochrones 
fitting of the upper MS indicate that Apsssw = 1-207 is sufficient for the Northeast Clump, 
but for R136 Apsssw = 1-525 is more appro priate. The same values have been obtained by 
applying a procedure similar to that used in |Pe Marchi et al.l ( 120111 ) 



It is commonly known that the steepness of the upper-MS and the paucity of massive 
stars makes it challenging to identify the MS turnoff (TO), and even harder to distinguish age 
differences of the order of few Myr, in star clusters younger than < 5 Myr using isochrone- 
fitting only. In such clusters a more robust age indicator is the luminosity o f the Turn-On 



(TOn), defined as the locus of the CMD where the PMS st ars join the MS f lStaufferill980 



Belikov et al.l Il998t iBaume et al.l 120031 ; ICignoni et al.l |2010|) . The age of a cluster can be 



defined as the time spent in the PMS phase by its most massive PMS star. 

The comparison between the CMDs of the two components shows that PMS stars in 
R136 are more than ~ 0.5 mag brighter compared to those found in the Northeast Clump. 
Since R136 is affected by a higher reddening than the Northeast clump, we conclude that 
PMS stars in R136 are intrinsically brighter, and therefore that in R136 SF lasted for at least 
1 Myr m ore than in the Northeas t clumrH It is possible to see the Northeast clump also in 
Fig. 11 of |Pe Marchi et al.l (120111 ). where the spatial distribution of PMS stars younger than 
4 Myr is compared to that of older objects with active mass accretion. 



2 We can exclude that the difference in magnitude between the two clumps is due to the distance, because 
AmFsssw — 1 would imply that the Northeast clump would be at ~ 79.5 kpc, far beyond even the Small 
Magellanic Cloud. 
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Discussion and Conclusions 



We used deep HST/WFC3 UV, optical and near-IR observations of the core of 30 Dor 
to study its stellar content. In all the filters we identified a dual structure that cannot be 
attributed to dust extinction. A careful analysis of the stellar content of the clumps indicates 
that, while the Northeast clump formed the majority of its stars between 2 and 5 Myr ago, 
R136 started to form stars likely ~ 2 Myr ago and was still active ~ 1 Myr ago. 

In retrospect, there are substantial indications in the literature that the Northeast 
Clump comprises a distinct, more evolved entity. Its two brightest member s, R141 and R142 



are early B s upergiants that must be older than the denizens of R136 itself (IFeast et al.lll960 



Parkerlll993l ). Many of the fainter late-0 /early- B stars, being evolved, are also sub ject to the 
same age argument ( Walborn &: Blades! 1 19971 ; iMassey &: Hunterlll998l ; iBosch et al.l 1999). In 
an analysis of star formation in the 30 Dor core based upon these observations, ISelmanl (119991 ) 
proposed three sequential events with mean ages less than 1.5 Myr, 2.5 Myr, and 5 Myr; 
their Figure 9 shows that the first is strongly concentrated toward R136, while the heaviest 
concentration of the intermediate and old objects coincides with the Northeast Clump. 

The age difference between R136 and the Northeast Clump can also explain the, until 
now, puzzling presence of two WC stars in 30 Dor, as these highly evolved objects are 
seldom found in giant Hll regions. One of them, M k 33d, is locat ed in an apparent multiple 
system at the inner edge of the Northeast clump (jMelnickill985l ). while the northward tail 
on the Northeast clum p encompasses the multiple system R140 containing the other WC 
along with two WN's (jMoffat et al.l 119871 ) . Finally, even the isolated M supergiant IR18 
([McGregor fc Hylandlll98ll ) somewhat further to the NE, long considered a field interloper, 
could actually be associated with the older clump. Positions of the mentioned stars with 
respect to the two sub-clusters are shown in Fig. HI 

This is, of course, a complex three-dimensional region, and some earlier spectral types 
also appear in the predo minantly older fields, partic ularly in the Mk 33 system. The 
02 If*/WN5 object Mk 30 JCrowther fc WalborrJl201ll ) that lies adjacent to R142 along the 
ridge of the Northeast clump could be associated with R136; but on the o ther hand, some 
such objects are also found in older young clusters (e.g., TS3 in NGC 2060; ISchild fc Testor 
19921 ). possibly as a result of binary evolution. In summary, the information presented here 
provides strong support for our completely independent recognition of a double system with 
an age offset between them in the core of 30 Dor. 

Characterizing the complex structure of massive star-forming regions such as NGC 2070 
is important to understand the early stages of cluster formation. The age difference found 
between R136 and the Northeast clump, as well as their morphology, may indicate that the 
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core of NGC 2070 is the result of a recent or ongoing merger between two su b-clusters. Inter- 



esting ly from the analysis of the radial velocities of apparently single O stars, iHenault-Brunet 



( 120121 . submitted) found evidence for an internal rotation in NGC 2070, that can be due to 
either a recent merger between the main core of NGC 2070 and a secondary cluster or by 
the agglomeration and dumpiness of the gas. 



The majority of the stars studied in iHenault-Brunetl (120121 ) coincides with R136, but 
5 sources coincide wit h the Northeast Clump. If the total mass of the system is ~ 10 5 M Q 
( lAndersen et al.ll2009l ) then the velocity difference, assuming a distance d = 4.5 pc, would 
be v — \J~(2GM/d) ~ 13Kms _1 for a hyperbolic encounter (i.e. energy), and ~ 9kms _1 



for a circular orbit. The maximum velocity difference reported by iHenault- Brunetl (120 121 ) 
is ~ 5kms _1 . If R136 and the Northeast Clump are undergoing a merger, then either we 
are observing the interaction with a low inclination (which would agree with the fact that 
we see the tidal tail on the plane of the sky), or the 2 systems are bound and are falling in 
with low velocity. A detailed analysis of the mass distribution in the two sub-clusters and 
a more sophisticated modeling of the two component interaction will be presented in two 
forthcoming papers (Sabbi et al. in prep, Gieles et al. in prep.) 

The presence of distinct components in a very y oung system is also consistent with recent 
studies of galactic giant molecular clouds (GMC) ( IHeyer et al.l 120090. which contradict the 
tradit ional paradigm that GMCs are gravitationally bound systems. iBallesteros-Paredes et al. 
( 1201 ll ) suggested that GMCs can undergo a hierarchical gravitational collapse, with the 
collapse occurring on scales ranging from individual cores to the whole cloud. In this sce- 
nario, because of the turbulent fragmentation of a GMC, SF will not be spre ad over the 



whol e GMC, but will be localized in gravitationally bound pockets of gas (1 Clark et al. 



20051 ; IClark. Bonnell. fc Klessenl I2008I ). Star clust ers will be the final products of the hi 



erarchical merger of such smaller sub-structures (jBonnell. Bate, fc Vine! l2003t iBatd 12009 
Federrath etaDboioh . 



We finally note that in a recent paper, iFujii et al.l (120121 ) suggested that the cuspy 
density profile of R136 and the large number of massive runaway stars that seem to escap e 
from this cluster are signatures of a post-core-collapse star cluster ( IMackey fc Gilmordl2003l ). 
Considering the small age of R136, such an early core collapse would occur only if the cluster 
had an ini tial density of p r > 10 6 M Q pc -3 , which is considerably higher than the current 



estimates. 



Fujii et al.l (120121 ) suggest that an efficient way to speed up the collapse process 



is to invoke the hierarchical merging of several smaller clusters. 

Observational evidence for a hierarchical process of SF co mes also from the complex 

struc ture of massive star-forming regions such as NG C 604 in M33 (iMafz-Apellaniz. Perez. &: Mas- Hesse 



20041 ). or NGC 346 in the Small Magellanic Cloud (ISabbi et al.l 120071 ). Also, iGennaro et al. 
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( 120 111 ) suggested that the elongation and mass segregation found in Westerlund 1 may be the 
footprint of a merger of multiple sub-clusters, formed almost coevally in the parental GMC. 
If massive star clusters form from mergers of smaller sub-systems, this can also explain the 
high fraction of rotating globular clusters found in our Galaxy A study of the radial veloci- 
ties and/or proper motions of the stars in the two clumps would provide a further test of the 
merging scenario. Such datasets are currently being obtained by us (e.g. Henault-Brunet et 
al. submitted; Lennon et al. in prep.). 



ES and DL are thankful to Chris Evans and Vincent Henault-Brunet for all the com- 
ments and suggestions that helped improve the letter. ES is also grateful to Aida Wofford and 
the Massive Stars and Starbursts Journal Club at STScI for the useful discussion and sugges- 
tions. The authors acknowledge support from grant GO-12499.01-A, and from NASA Hubble 
Fellowship grant HST-HF-51270.01-A to SdM, awarded by STScI, operated by AURA for 
NASA, contract NAS 5-26555. N.P. acknowledges partial support by STScI-DDRF grant 
D0001.82435. 



REFERENCES 

Andersen, M., Zinnecker, H., Moneti, A., McCaughrean, M.J., Brandl, B., Brandner, W., 
Meylan, G., Hunter, D. 2009, ApJ, 707, 1347 

Anderson, J., Bedin, L.R., Piotto, G., Yadav, R. S., & Bellini, A. 2006, A&A, 454, 1029 

Anderson, J., & King, I.R. 2006, PSFs, Photometry and Astrometry for the ACS/WFC 
(ISR-ACS 2006-01; Baltimore: STScI) 

Ballesteros-Paredes, J. Hartmann, L.W., Vazquez-Semadeni, E., Heitsch, F., Zamora-Aviles, 
M.A. 2011, MNRAS, 277, 362 

Bastian, N. Goodwin, S.P 2006, MNRAS, 369, L9 

Bate, M.R. 2009, MNRAS, 392, 590 

Baume, G., Vazquez, R.A., Carraro, G., Feinstein, A. 2003, A&A, 402, 549 

Belikov, A.N., Hirte, S., Meusinger, H., Piskunov, A.E., Schilbach, E. 1998, A&A, 332, 575 

Bellini, A., Anderson, J., & Bedin, L.R. 2011, PASP, 123, 622 

Bellini, A., & Bedin, L.R. 2009, PASP, 121, 1419 



-8 - 

Bonnell, LA., Bate, M.R., Vine, S.G. 2003, MNRAS, 343, 413 
Bosch, G., Terlevich, R., Melnick, J., Selman, F. 1999, A&AS, 137, 21 
Brandl, B.R. et al. 2004, ApJS, 154, 188 

Campbell, M.A., Evans, C.J., Mackey, A.D., Gieles, M., Alves, J., Ascenso, J., Bastian, N., 
Longmore, A.J. 2010, MNRAS, 405, 421 

Clark, P.C, Bonnell, I.A., Zinnecker, H., Bate, M.R. 2005, MNRAS, 359, 809 

Clark, P.C, Bonnell, LA. & Klessen, R.S. 2008, MNRAS, 386, 3 

Cignoni, M., Tosi, M., Sabbi, E., Nota, A., Degl'Innocenti, S., Prada Moroni, P.G., Gal- 
lagher, J.S. 2010, 712, L63 

Crowther, P.A., et al. 2010, MNRAS, 405, 421 

Crowther P.A., & Walborn, N.R. 2011, MNRAS, 416, 1311 

de Kother, A., Heap, S.R, & Hubeny, I. 1998, ApJ, 509, 879 

De Marchi, G. et al. 2011, ApJ, 739, 27 

Elson, R.A., Fall, S.M., Freeman, K.C. 2987, ApJ, 323, 54 

Feast, M.W., Thackeray, A.D., & Wesselink, A.J. 1960, MNRAS, 121, 337 

Federrath, C, Banerjee, R., Clark, P.C, Klessen, R.S. 2010, ApJ, 713, 269 

Fujii, M.S., Saitoh, T.R., & Portegies Zwart, S.F. 2012, ApJ, 753, 85 

Gennaro, M., Bradner, W., Stolte, A, Hemmimg, Th. 2011, MNRAS, 412, 2469 

Heckman, T.M., Kauffmann, G., Brinchmann, J., Chariot, S., Tremonti, C, White, S.D.M. 
2004, ApJ, 613, 109 

Helling, C, Winters, J.M., Sedlmayr, E. 2000, A&A, 358, 651 

Henault-Brunet, V., Gieles, M., Evans, C.J., Sana, H., Bastian, N., Maiz-Apellaniz, J., Tay- 
lor, W.D., Markova, N., Bressert, E., de Koter, A., van Loon, J. Th. 2012, submitted. 

Heyer, M., Krawczyk, C, Duval, J., Jackson, J.M. 2009, ApJ 699, 1092 

Kennicutt, R.C. Jr. 1984, ApJ, 287, 116 



-9 - 

Mackey, A.D, & Gilmore, G.F. 2003, MNRAS, 338, 85 
Mafz-Apellaniz, J., Perez, E., Mas-Hesse, J.M. 2004, AJ, 128,1196 

Marigo. P., Girardi, L., Bressan, A., Groenwegen, M.A.T., Silva, L., Granato, G.L. 2008, 
A&A 482, 883) 

Massey, P., & Hunter, D.A. 1998, ApJ, 493, 180 

McGregor, P. J., & Hyland, A.R. 1981, ApJ, 250, 116 

Melnick, J. 1985, A&A, 153, 235 

Meuer, G.R., Heckman, T.M., Lehner, M.D., Leitherer, C, Lowenthal, J. 1997, AJ, 114, 54 

Meylan, G. 1993, in Smit G.H., Brodie J.P, eds., ASP Conf. Srr. Vol. 48, The Globular 
Cluster-Galaxy Connection. Astron. Soc. Pac, San Francisco, p. 588 

Moffat, A.F.J., Niemela, V.S., Phillips, M.M., Chu, Y.-H., & Seggewiss, W. 1987, ApJ, 312, 
612 

Nikolaev S., Drake, A.J., Keller, S.C., Cook, K.H., Dalai, N., Griest, K., Welch, D.L., 
Kanbur, S.M. 2004, ApJ, 601, 260 

Parker, J.W. 1993, AJ, 106, 560 

Sabbi, E., Sirianni, M., Nota, A., Tosi, M., Gallagher, J.S, Meixner, M., Oey, M.S., Walter- 
bos, R., Pasquali, A., Smith, L.J., Angeretti, L. 2007, AJ, 133, 2430 

Schaefer G.H., Simon, M., Prato, L., Barman, T. 2008, AJ, 135, 1659 

Schild, H., & Testor, G. 1992, A&AS, 92, 729 

Selman, F., Melnick, J., Bosch, G., Terlevich, R. 1999, A&A, 347, 532 
Shapley, A.E., Steidel, C.C., Pettini, M., Adelberger, K.L. 2003, ApJ, 588, 65 
Stauffer, J.R. 1980, AJ, 85, 1341 

Tognelli, E., Prada Moroni, P.G., Degl'Innocenti, S. 2011, A&A, 533A, 109 
Walborn, N.R., & Blades, J.C. 1997, ApJS, 112, 457 



This preprint was prepared with the A AS IATgX macros v5.2. 



- 10 - 




Fig. 1.— CMDs of ^F336W VS- ?TlF336W — ""T-F438W (panel A), mF438W V S- ?^F438W — ?™F555W 

(panel B), and mps^w vs - m F555W — ^fsuw (panel C), of the stars found in the NGC 2070 re 



gion with photometric errors < 0.15 mag in both filters. Pisa isochrones (solid; iTognelli et al. 
20 111 ) for PMS stars of metallicity Z=0.008 and ages=l, 2 and 5 Myr are sho wn in blue 



green and red respectively. The dashed isochrones are from iMarigo et al. I (120081 ). A 10 Myr 
old stellar population with Z=0.008 is shown in cyan. A 4 Gyr old stellar population with 
Z= 0.004 and two distinct values of reddening (E(B-V)=0.14 and 0.4) is shown in magenta 
and dark green respectively. 
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Fig. 2. — Stellar density charts of the stars found in the core of NGC 2070 in the filters 
F336W , F438W, F555W, F814W and F160W. In each map the highest level is 5 times 
higher than the lower one. The lower-right plot shows that ratio between the F160W and 
the F555W maps. Isodensity contours from the F555W map are plotted for reference. In all 
the maps North is up, and East is to the left. 




Fig. 3. — UVIS image of 30 Dor in the F555W filter. North is up, East to the left. Isodensity 
contours as derived in Fig. [2] are shown for reference. 




Fig. 4. — F555W UVIS image zoomed on the center of 30 Dor. The dashed circle and 
rectangle indicate the regions used to derive the CMDs shown in Fig. [51 Stars mentioned in 
the Discussion are highlighted. Two isodensity contours are shown to guide the eye. 
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Fig. 5. — CMDs of m.F555w vs. mpsssw — ^F8i4W of the stars with photometric error < 0.15 
mag in both the F555W and F814W filters found within ~ 10" from the center of R136 
(left panel) and in a ~ 14.4" x 27" box centered on the Northeast clump (right panel). Pisa 
isochrones for PMS stars with metallicity Z=0.008 and an age of 1, 2 and 5 Myr are plotted 
in blue, red and green respectively. In both the CMDs a distance modulus of 18.5 was 
assumed. R136 is better reproduced assuming an Ap555w = 1.525. For the Northeast clump 
A-F555W — 1-207 is sufficient. 



